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Abstract 

The Schrodinger equation with the hyperbolic Scarf potential reported so far in the literature is 
somewhat artificially manipulated into the form of the Jacobi equation with an imaginary argument 
and parameters that are complex conjugate to each other. Instead we here solve the former equation 
anew and make the case that it reduces straightforward to a particular form of the generalized real 
hypergeometric equation whose solutions are referred in the mathematics literature as the finite 
Romanovski polynomials in reference to the observation that for any parameter set only a finite 
number of such polynomials appear orthogonal. This is a qualitatively new integral property that 
does not copy none of the features of the Jacobi polynomials. In this manner the finite number of 
bound states within the hyperbolic Scarf potential is brought in correspondence to a finite system 
of orthogonal polynomials of a new class. 

This work adds a new example to the circle of the problems on the Schrodinger equation. The 
techniques used by us extend the teachings on the Sturm-Liouville theory of ordinary differential 
equations beyond their standard presentation in the textbooks on mathematical methods in physics. 

La solution a la ecuacion de Schrodinger con el potential de Scarf hiperbolico reportada hasta ahora 
en la literatura fisica esta manipulada artificialmente para obtenerla en la forma de los polinomios 
de Jacobi con argumentos imaginarios y parametros que son complejos conjugados entre ellos. En 
lugar de eso, nosotros resolvimos la nueva ecuacion obtenida y desarrollamos el caso en el que 
realmente se reduce a una forma particular de la ecuacion hipergeometrica generalizada real, cuyas 
soluciones se refieren en la literatura matemtica como los polinomios finitos de Romanovski. La 
notation de finito se refiere a que, para cualquier parametro fijo, solo un mimero finito de dichos 
polinomios son ortogonales. Esta es una nueva propiedad cualitativa de la integral que no surge 
como copia de ninguna de las caractersticas de los polinomios de Jacobi. De esta manera, el mimero 
finito de estados en el potential de Scarf hiperbolico es consistente en correspondencia a un sistema 
finito de polinomios ortogonales de una nueva clase. 

PACS numbers: 02.30.Gp, 03.65. Ge, 12.60.Jv 
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I. INTRODUCTION 



The exactly solvable Schrodinger equations occupy a pole position in quantum mechanics 
in so far as most of them relate directly to physical systems. Suffices to mention as prominent 
examples the quantum Kepler-, or, Coulomb problem and its importance in the description 
of the discrete spectrum of the hydrogen atom tlfl, the harmonic-oscillator, the Hulthen, 

n n 

and the Morse potentials with their relevance to vibrational spectra [2|, |3|. Another good 
example is given by the Poschl- Teller potential 41] which appears as an effective mean field in 
many-body systems with 5-interactions [5( . In terms of path integrals, the criteria for exact 
resolvability of the Schrodinger equation is the existence of exactly solvable corresponding 
path integrals jf]]. 

There are various methods of finding the exact solutions of the Schrodinger equation (SE) 
for the bound states, an issue on which we shall focus in the present work. The traditional 
method, to be pursued by us here, consists in reducing SE by an appropriate change of the 
variables to that very form of the generalized hypergeometric equation [7| whose solutions 
are polynomials, the majority of them being well known. The second method suggests to 
first unveil the dynamical symmetry of the potential problem and then employ the relevant 
group algebra in order to construct the solutions as the group representation spaces 0, S]. 
Finally, there is also the most recent and powerful method of the super-symmetric quantum 
mechanics (SUSYQM) which considers the special class of Schrodinger equations (in units 

nn 

of h = 1 = 2m) that allow for a factorization according to [10(-[12j, 

H(z) = A + (z)A~(z) + e , 

A±{z) = (±4z + U(z)) • (1) 



dz 

Here, H(z) stands for the (one-dimensional) Hamiltonian, U(z) is the so called super- 
potential, and A ± (z) are the ladder operators connecting neighboring solutions. The super- 
potential allows to recover the ground state wave function, ip gst (z), as 

iP gst (z)~e-f Zu ^ d y. (2) 

The excited states are then built up on top of ipgst{z) through the repeated action of the 
A + (z) operators. 



A. The trigonometric Scarf potential. 



The super-symmetric quantum mechanics manages a family of exactly solvable potentials 
(see Refs. [ll| -[l3| for details) one of which is the so called trigonometric Scarf potential 
(Scarf I) 14j . here denoted by v t (z) and given by 



-a 2 + (a 2 + b 2 — aa) sec 2 az — b{2a + a) tan az sec az . 



(3) 



It has been used in the construction of a periodic potential and employed in one-dimensional 
crystal models in solid state physics. 

The exact solution of the Schro ding er equation with the trigonometric Scarf potential 



(displayed in Fig. 1) is well known 

P£> a (x) : as 
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13] and given in terms of the Jacobi polynomials, 



ip n (x 
w 1 - 



y/(l-x)7(l + x) S P2 ^ '(z), 

2{ X ) = (1 - xY'^il + x) s ~ 



X 
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smaz 
1 



a 



a — b), 



a 



{a + b) 



(4) 



x) stands for the weight function from which the Jacobi polynomials 



Here, w 1 2> u 
' ' 9 (x) are obtained via the Rodrigues formula. 
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Fig. 1. The trigonometric Scarf potential (Scarf I) for the set of parameters, a = 10, b = 5, and 
a = 1. The horizontal lines represent the discrete levels. 



The corresponding energy spectrum is obtained as 



e n = e n + a 2 = (a + na) 2 



(5) 



B. The hyperbolic Scarf potential. 

By means of the substitutions 



-ia , b — > b , 



(6) 



Scarf I is transformed into the so called hyperbolic Scarf potential (Scarf II), here denoted 
by v^' b \z) and displayed in Fig. 2, 



(0,6)/ \ 
v h ( z ) 



ar + 



a 2 — aa) sech 2 az + b(2a + a) sech azt&nhaz . 



(7) 



The latter potential has also been found independently within the framework of super- 
symmetric quantum mechanics while exploring the super-potential 

U(z) = a tanh az + b sech az . (8) 

Upon the above substitutions and in taking a = 1 for simplicity the energy changes to 

£n = e n — a 2 = — (a — n) 2 , n — 0, 1, 2, ... < a . (9) 



V (I) 




Fig. 2. The hyperbolic Scarf potential (Scarf II) for the set of parameters, a = 10, 6 = 5, and 
a = 1. The horizontal lines represent the energies, e n , of the bound states. 

It is important to notice that while the trigonometric Scarf potential allows for an infinite 
number of bound states, the number of discrete levels within the hyperbolic Scarf potential is 
finite, a difference that will be explained in section III below. Yet, the most violent changes 
seem to be suffered by the Jacobi weight function in eq. (BJ and are due to the opening of 
the finite interval [—1, +1] toward infinity, 

x = sinaz G [—1, 1] — ► x = sinhaz G [—00, +00] . (10) 



In this case, the wave functions become 



$ n (-ix) = (l + x 2 )- a *e- bt ™- lx c n Pf"(-ix), v = ib-a- 1 -. (11) 

Here, c n is some state dependent complex phase to be fixed later on. The latter equation 
gives the impression that the exact solutions of the hyperbolic Scarf potential rely exclusively 
upon those peculiar Jacobi polynomials with imaginary arguments and complex indices. We 
here draw attention to the fact that this needs not be so. 
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C. The goal. 



The goal of this work is to solve the Schrodinger equation with the hyperbolic 
Scarf potential anew and to make the case that it reduces in a straightforward 
manner to a particular form of the generalized real hypergeometric equation 
whose solutions are given by a finite set of real orthogonal polynomials. In this 
manner, the finite number of bound states within the hyperbolic Scarf potential 
is brought in correspondence with a finite system of orthogonal polynomials of 
a new class. 

These polynomials have been discovered in 1884 by the English mathematician Sir Edward 
John Routh 18] and rediscovered 45 years later by the Russian mathematician Vsevolod 
Ivanovich Romanovski in 1929 [l9| within the context of probability distributions. Though 
they have been studied on few occasions in the current mathematical literature where they 

nn 

are termed to as "finite Romanovski" [20|-[23], or, "Romanovski-Pseudo-Jacobi" polyno- 



mials 24J, they have been completely ignored by the textbooks on mathematical methods 



in p hysics, and surprisingly enough, by the standard mathematics textbooks as well 



a, 



25j-[28|. The notion "finite" refers to the observation that for any given set of parameters 
(i.e. in any potential), only a finite number of polynomials appear orthogonal. 

The Romanovski polynomials happen to be equal (up to a phase factor) to Jacobi polyno- 
mials with imaginary arguments and parameters that are complex conjugate to each other, 
much like the sinhz = isiniz relationship. Although one may (but has not to) deduce the 
local characteristics of the latter such as generating function and recurrence relations from 
those of the former, the finite orthogonality theorem is qualitatively new. It does not copy 
none of the properties of the Jacobi polynomials but requires an independent proof. 

Our work adds a new example to the circle of the typical quantum mechanical problems 



291 ] . The techniques used by us here extend the study of the Sturm-Liouville theory of 
ordinary differential equations beyond that of the usual textbooks. 

A final comment on the importance of the potential in eq. ©. The hyperbolic Scarf 
potential finds various applications in physics ranging from electrodynamics and solid state 
physics to particle theory. In solid state ph ysic s Scarf II is used in the construction of 
more realistic periodic potentials in crystals [30] than those built from the trigonometric 
Scarf potential. In electrodynamics Scarf II appears in a class of problems with non-central 

6 



potentials (see section IV ). In particle physics Scarf II finds application in studies of the non- 
perturbative sector of gauge theories by means of toy models such as the scalar field theory 
in (1+1) space-time dimensions. Here, one encounters the so called "kink -like" solutions 
which are no more but the static solitons. The spatial derivative of the kink-like solution 
is viewed as the ground state wave function of an appropriately constructed Schrodinger 
equation which is then employed in the calculation of the quantum corrections to first order. 
In Ref. [l3| it was shown that specifically Scarf II is amenable to a stable renormalizable 
scalar field theory. 

The paper is organized as follows. In the next section we first highlight in brief the basics 
of the generalized hypergeometric equation and then relate it to the Schrodinger equation 
with the hyperbolic Scarf potential. The solutions are obtained in terms of finite Romanovski 
polynomials and are presented in section III. Section IV is devoted to the disguise of the 
Romanovski polynomials as non-spherical angular functions. The paper ends with a brief 
summary. 

II. MASTER FORMULAS FOR THE POLYNOMIAL SOLUTIONS TO THE GEN- 
ERALIZED HYPERGEOMETRIC EQUATION. 

All classical orthogonal polynomials appear as solutions of the so called generalized hy- 



pergeometric equation (the presentation in this section closely follows Ref. [221 ]) 



a{x)y';{x)+r{x)y' n {x)-X n y n {x) = 0, (12) 
cr(x) = ax 2 + bx + c, t(x) = xd + e , A n = n{n — l)a + nd . (13) 



There are various methods for finding the solution, here denoted by 

y n (x) = P n 

I A p ^ 

with the symbol P n 



Id e 






X 


\ a b c 





(14) 



x j in which the equation parameters have been made explicit 
standing for a polynomial of degree n, A„ being the eigenvalue parameter, and n = 0, 1, 2, .... 



Id e 






X 


\ a b c 





In Ref. 



22J a master formula for the (monic, P n ), polynomial solutions has been derived by 
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Koepf and Masjed-Jamei, according to them one finds 

d e 
a b c 



Pr, 



G^\a,b, c,d,e)x k , 



G 



(n) 




(k - n), 



2ae—bd i -i d 

2aVb 2 -4ac ^ 2a 



n 



2n 



2^b 2 - 4ac 



b + \Jb 2 - Aac 



(15) 



Though the formal proof of this relation is bit lengthy, its verification with symbolic math- 
ematical softwares like Maple is straightforward. The a = case is treated as the a — > 
limit of eq. (115j) and leads to 2-P0 i n place of 2-^1- Notice that G% are n °t normalized. On 
the other side, eq. (1121) can be treated alternatively as described in the textbook by Nikiforov 
and Uvarov in Ref. Q where it is cast into a self-adjoint form and its weight function, w(x), 
satisfies the so called Pearson differential equation, 

d 



dx 



(a(x)w(x)) = t(x)w(x) 



(16) 



The Pearson equation is solved by 

W 



w[x) 



d e 
a b c 



x 



exp 



(d - 2a)x + {e-b) 



bx + c 



(17) 



and shows how one can calculate any weight function associated with any parameter set of 
interest (we again used a symbol for the weight function that makes explicit the parameters 
of the equation). The corresponding polynomials are now classified according to the weight 
function and are built up from the Rodrigues representation as 



Pn 



Id e 






X 


\ a b c 





Tl k k zUd +(n + k- 2)a)P n 



Id e 






X 


\ a b c 





1 



w 



Id e 






X 


\ a b c 





X 



d n 
dx n 



{ax 2 + bx + c) n W 



Id e 


\ 

X 


') 


\ a b c 







(18) 



The master formulas in the respective eqs. (11 5p . and ([TBI allow for the construction of all the 
polynomial solutions to the generalized hypergeometric equation. One identifies as special 
cases the canonical parameterizations known as 
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the Jacobi polynomials with a = — 1, b — 0, c = 1, d — —7 — 5 — 2, and e = —7 + 5, 



• the Laguerre polynomials with a = 0, b = 1, c = 0, d = — 1, and e = a + 1, 

• the Hermite polynomials with a = b = 0,c = l,d = —2, and e = 0, 

• the Romanovski polynomials with a — 1, b = 0, c = 1, d = 2(1 — p), and e = q with 

p> 0, 

• the Bessel polynomials with a — 1, 6 = 0, c = 0, d — a + 2, and e = (3. 

All other parameterizations can be reduced to one of the above five by an appropriate shift of 
the variables. The first three polynomials are the only ones that are traditionally presented 
in the standard textbooks on mathematical methods in physics such like 25j- 28), while the 
fourth and fifth seem to have escaped due attention. Notice, the Legendre, Gegenbauer, and 
Chebychev polynomials appear as particular cases of the Jacobi polynomials. The Bessel 
polynomials are not orthogonal in the conventional sense, i.e. within a real interval, and 
will be left out of consideration. 

Some of the properties of the fourth polynomials have been studied in the specialized 
mathematics literature such as Refs. [2^], [21], 23]. Their weight function is calculated from 
eq. (TI7I) as 

w M(x) = (x 2 + i)-P e i^-^. (19) 
This weight function has first been reported by Routh 18], and independently Romanovski 



[l9l | . The polynomials associated with eq. ([191) are named after Romanovski and will be 
denoted by Rm' q \x). They have non-trivial orthogonality properties over the infinite interval 
[—00, +00]. Indeed, as long as the weight function decreases as x~ 2p , hence integrals of the 
type 

/+00 
w {p ' q \x)R^\x)R^} q \x)dx (20) 
-00 

are convergent only if 

m + m'<2p-l, (21) 



meaning that only a finite number of Romanovski polynomials are orthogona' 
the reason for the term "finite" Romanovski polynomials (details are given Ref. 
differential equation satisfied by the Romanovski polynomials reads as 



This is 



3l|]). The 



2 p(p,<?), 



d 2 R 



x 



d 2 x 



+ (2(-p + l)x + q) 



dRt q) ( 
dx 
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(n(n-l)+2n(l-p))R { ™\x) = 0. (22) 



In the next section we shall show that the Schrodinger equation with the hyperbolic Scarf 
potential reduces precisely to that very eq. (1221) . 



A. The real polynomial equation associated with the hyperbolic Scarf potential. 

The Schrodinger equation for the potential of interest when rewritten in a new variable, 



33(, taken by us as 



x, introduced via an appropriate point canonical transformation 
z = f(x) = sinh -1 x, is obtained as: 

(1 +x ^fiM + Jj£) + (Z^±i) _ K^±l) x + £ „) 9 „ W = o , ( 23 ) 

with g n (x) = ?/v(sinh -1 x), and e n = e n — a 2 . Equation f|T9|) suggests the following substitu- 
tion in eq. f l23l ) 

gjx) = (l + x 2 )2 e ~^ tan ' lx D^' a \x), x = sinhz, -oo < x < +00. (24) 
In effect, eq. (|23l) reduces to the following equation for Dn (x), 

(1 + J& +((2{i + l)x-a) 



dx 2 dx 



(25) 

If the potential equation (|25|) is to coincide with the Romanovski equation (122]) then 

• first the coefficient in front of the l/(x 2 + 1) term in (125j) has to vanish, 

• the coefficients in front of the first derivatives have to be equal, i.e. 2(— p + 1) + q — 
(2(3 + l)x-a, 

• the eigenvalue constants should be equal too, i.e. e n + j3 2 = —n((n — 1) + 2(1 — p)). 

The first condition restricts the parameters of the D^' a \x) polynomials to 

a 2 

a + a 2 -b 2 + ^ + (3-l3 2 = 0, (26) 
-b-2ab+--a(3 = 0. (27) 
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Solving the equations (J2BI) . f[2"T|) . for a and /3 results in 

@ = -a, a = 26. (28) 

The second condition relates the parameters a and j3 to p, and q, and amounts to 

(3 = -a = -p+^, - a = q = -2b. (29) 

Finally, the third restriction leads to a condition that fixes the Scarf II energy spectrum as 

e n = -(a-n) 2 . (30) 

In this way, the polynomials that enter the solution of the Schrodinger equation will be 

D^=-^ b \x) = Rt a+hq= - 2b \x). (31) 

They are obtained by means of the Rodrigues formula from the weight function w^ a+ ^'^ 2b \x) 
as 

1 1 rt m 

(a+ 5 ,-26), v _ _fl I r 2\m (o+i,-2&)/ \ 

^ ^ " w (a + I,-26) (;E)c / a;m ^ + X ^ W' 

As a result, the wave function of the nth level, ip n , takes the form 

^ n (^ = sinh- 1 a;) := g n {x) = 1 _ \/ {1 + x 2 y( a+ ^ e- 2 ^ 1 -R^'^jx), 

/ d sinh 1 x 
V <ix 

(isinh -1 x = dx . (33) 

VTTx^ 

The orthogonality integral of the Schrodinger wave functions gives rise to the following 
orthogonality integral of the Romanovski polynomials, 

"+oo /*+oo 1 . 

rh(z)M*)dz= / (l + ^-^^e- 26 ^ 1 ^^ 5 '^^)^ 5 '^^^)^, (34) 



CO 



which coincides in form with the integral in eq. ( 1201) and is convergent for n < a. That only 
a finite number of Romanovski polynomials are orthogonal, is reflected by the finite number 
of bound states within the potential of interest, a number that depends on the potential 
parameters in accord with eq. (I2T]) . 

As to the complete Scarf II spectrum, it has been constructed in Ref. [9J within the 
dynamical symmetry approach [8]. There, the Scarf II potential has been found to possess 
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SU(1, 1) as a symmetry group. The bound states have been assigned to the discrete unitary 
irreducible representations of SU(1, 1). The scattering and resonant states (they are beyond 
the scope of the present study) have been related to the continuous unitary and the non- 
unitary representations of 577(1, 1), respectively 

A comment is in place on the relation between the Romanovski polynomials and the 
Jacobi polynomials of imaginary arguments and parameters that are complex conjugate to 
each other. Recall the real Jacobi equation, 

fj2pj,S( T \ dP~*> 5 (r) 
(1 ~ %2) dx 2 + (7 " 6 ~ (7 + 6 + 2)X) ^x -< n + ^ + 5 + l ) P n ( x ) = • (35) 

As mentioned above, the real Jacobi polynomials are orthogonal within the [—1, 1] interval 
with respect to the weight-function in eq. (HI). Transforming to complex argument, x — > ix, 
and parameters, 7 = 5* = c + id, eq. ( l35i) transforms into 

j2 pc+id,c—idf- \ J pc+id,c-id( ■ \ 

(1 + x 2 ) n J - [tX) + (-2d + 2(c+l)x) a " 



dx 2 dx 

+ n(n + 2c + l)P^ +id ' c - id (ix) = . (36) 

Correspondingly, the weight function turns to be 

w c+id ' c - id {ix) = (1 + x 2 ) c e- Mtan_1 x , (37) 

and it coincides with the weight function of the Romanovski polynomials in eq. (fT9l) upon 
the identifications c = —p, and q = —2d. This means that p^ ld ^- ld [ix) will differ from the 



Romanovski polynomials by a phase factor found as i n in Ref . 35( , among others 



i n Pn 



2(1 — p) iq 
-1 1 



ix 




x 



(38) 



Because of this relationship the Romanovski polynomials have been termed to as 
"Romanovski-Pseudo- Jacobi" by Lesky [24j. The relationship in eq. (|38|) tells that the 
Rn' 9 \x) properties translate into those of P~ p ~ l i~ v+l i(ix) and visa versa, and that it is a 
matter of convenience to prefer the one polynomial over the other. When it comes up to 
recurrence relations, generating functions etc. it is perhaps more convenient to favor the 
Jacobi polynomials, though the generating function of the Romanovski polynomials is also 
equally well calculated directly from the corresponding Taylor series expansion [3jj. How- 
ever, concerning the orthogonality integrals, the advantage is clearly on the side of the real 
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Romanovski polynomials. This is so because the complex Jacobi polynomials are known 
for their highly non-trivial orthogonality properties which depend on the interplay between 



integration contour and parameter values 36j. For this reason, in random matrix theory 



371 ] the problem on the gap probabilities in the spectrum of the circular Jacobi ensemble 
is treated in terms of the Cauchy random ensemble, a venue that heads one again to the 
Romanovski polynomials (notice that for p = 1, q = the weight function in eq. (1191) reduces 
to the Cauchy distribution). 

In summary, and for all the reasons given above, the Romanovski polynomials qualify 
as the most adequate real degrees of freedom in the mathematics of the hyperbolic Scarf 
potential. 

III. THE POLYNOMIAL CONSTRUCTION 

The construction of the i?i + 2 ' (x) polynomials is now straightforward and based upon 
the Rodrigues representation in eq. (ITS]) where we plug in the weight function from eq. (jl~9l) . 
In carrying out the differentiations we find the lowest four (unnormalized) polynomials as 

R\ 2 ' \x) = 1, (39) 
R [ a+1 ^- 2b ) ( x ) = _2b + (1 - 2a)x , (40) 
R i a+ h- 2b ) ( x ) = 3 _ 2a + 452 _ g^j _ a } x + ( 6 _ 10fl + 4a 2) x 2 ^ ^ 

R { ^'~ 2b \x) = -266 + Ylah - 8b 3 + [-3(-15 + 16a - 4a 2 ) + 12(3 - 2a)6 2 ]x 

+ (-726 + 84a6 - 24a 2 6)x 2 + 2(-2 + a) (-15 + 16a - 4a 2 )x 3 . (42) 

As illustration, in Fig. 3 we show the Scarf II wave functions of the first and third levels. 



The finite orthogonality of the Romanovski polynomials becomes especially transparent 
in the interesting limiting case of the sech 2 z potential (it appears in the non-relativistic 
reduction of the sine- Gordon equation) where one easily finds that the normalization con- 
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Fig. 3. Wave functions for the first and third levels within the hyperbolic Scarf potential, 
stants, Ni a+2 '°\ are given by the following expressions: 

/ M,)y = H.-1W-D a> 

(n^X = ^-i)r(°-2) _ a > 
V / r («-|) 

[nI 2 ' ) = v y — J -^- ! -(Aa 2 - 16a + 15) 2 , a > 3 etc. (43) 

Software like Maple or Mathematica are quite useful for the graphical study of these func- 
tions. The latter expressions show that for positive integer values of the a parameter, a = n, 
only the first (n — 1) Romanovski polynomials are orthogonal (the convergence of the inte- 
grals requires n < a), as it should be in accord with eqs. ( 12T1) . ([91). The general expressions 
for the normalization constants of any Romanovski polynomial are defined by integrals of 
the type f_ °°(1 + x 2 )~ ( - p ~ n ' , e qtan lx dx and are analytic for (p — n) integer or semi-integer. 

IV. ROMANOVSKI POLYNOMIALS AND NON-SPHERICAL ANGULAR 
FUNCTIONS IN ELECTRODYNAMICS WITH NON-CENTRAL POTENTIALS. 



In recent years there have been several studies of the bound states of an electron within a 



compound Coulomb- and a non-central potential (see Refs. 38 
Let us assume the following potential 

V 2 (6) 



391 ] and references therein). 



V(r,6)=V c (r) + 



V 2 (9) 



c cot 9 . 



(44) 
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where Vc(r) denotes the Coulomb potential and 9 is the polar angle (see Fig. 4). The 
corresponding Schrodinger equation reads 



1 d o d 



1 d d 
smu- 



1 



2 



r 2 dr dr r 2 sin 9 89 89 r 2 sin 2 9 dip 2 
and is solved as usual by separating variables, 



V(r,6) 



*{r,0,<p) = E*(r,e,<p), 

(45) 



(46) 



The radial and angular differential equations for TZ(r) and 0(9) are then found as 
d 2 K{r) 2 dK(r) 



dr 2 



r dr 



(Vc(r)+E)- l «±± 



and 



d 2 9(^) , n ,de(6) 



d9 2 



d6 



1(1 + 1) - V 2 (9) - 



m 



sin^ 



TZ(r) = 0, 



0(9) =0. 



(47) 



(48) 



with 1(1 + 1) being the separation constant. From now on we will focus attention on eq. f|48l) . 
It is obvious that for V 2 (9) = 0, and upon changing variables from 9 to cos 9, it transforms 
into the associated Legendre equation. Correspondingly, 6(6*) approaches the associated 
Legendre functions, 

9(0) V2 ^° P™(cos8) , (49) 
an observation that will become important below. 




Fig. 4. The non-central potential V(r,6), here displayed in its intersection with the x = 
plane, i.e. for r = \/y 2 + z 2 , and 6 = tan -1 |. The polar angle part of its exact solutions is 
expressed in terms of the Romanovski polynomials. 
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In order to solve eq. (|48p we follow the prescription given in 38j and begin by substituting 
the polar angle by the new variable, z, introduced via 9 — > f(z). This transformation leads 
to the new equation 



dz 2 



+ 



/"(*) 
/'(*) 



+ f(z) cot f(z) 



dz 



m 



+ -r 2 (/(-)) + /(/ + i)- , ;, „, , /'%:) ij(z) =o, 

(50) 

with f'(z) = and if>(z) defined as ip{z) = Q(f(z)). Next, one requires the coefficient 

in front of the first derivative to vanish which transforms eq. (1501) into a Id Schrodinger 
equation. This restricts f(z) to satisfy the differential equation 

f"(z) 



f(z) oot/(z) 



(51) 



which is solved by f(z) = 2 tan 1 e z . With this relation and after some algebraic manipula- 
tions one finds that 



sin 9 



1 



cosh z ' 



cos 9 



tanh^ 



(52) 



and consequently, 

f'(z) = sin f(z) = sechz. (53) 

Equation ( |52l) implies sinh^; = — cot^, or, equivalently, 9 = cot _1 (— sinh z). Upon substi- 
tuting eq. ( |53l) into eqs. (jHj), and ( 15D1 . one arrives at 



d 2 j){z) 
dz 2 



+ 



1(1 + 1) 



1 , 1 

+ ctanhz 



cosh z 



cosh z 



m 



ip(z) := 0(9 = cot' 1 (- sinh z)) , 0(9) := ip (z = sinh -1 (- cot 9)) 



clef 



4>(z) = 0, 
(54) 



Taking in consideration eqs. f!7l). (H4]) . and (152^) one realizes that the latter equation is pre- 
cisely the one-dimensional Schrodinger equation with the hyperbolic Scarf potential where 



1(1 + 1) = -(b 2 - a (a + 1)) , c=-b(2a + l), 



m 



(a — n)' 



m > . 



(55) 



The two parameters of the Romanovski polynomials have to be determined from the system 
of the last three equations, meaning that the I, m, and c constants can not be independent. 
There exist various choices for a and b. If defined on the basis of the first two equations, 
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one encounters 




(56) 

Substitution of a into the third equation imposes a constraint on I as a function of m, c, 
and n. A second choice for a and b is obtained by expressing a from the third equation in 
terms of m, and n as a = m + n and substituting in the second equation to obtain b as 

c 



2(m + n) + 1 

Then the first equation imposes the following restriction on I 



(57) 



X-=(b 2 -a(a+l)), l = -l + j± + X. (58) 

This I value which is not necessarily integer, is the one that enters the well known energy, 
E nr i = —Z 2 e 2 fi/ (2fr 2 (n r + / + l) 2 ), attached to the radial solution, thus leading to a (discrete) 
spectrum that no longer bears any resemblance to the 0(4) degeneracy. This is the path 
pursued by Ref. [381 ] . We here instead take a third chance and express a, b, and c as functions 
of I alone according to 

a = b = 1(1 + 1), n = a - m = 1(1 + 1) - m, c= -b(2a + l). (59) 

This choice allows to consider integer / values. 

In making use of Eqs. (|2"9l . (T5TT) . the solution for 9 becomes 

6(0) = tpn=i(i+i)-m (z = sinh _1 (-cot^)) 

= (1 + cot 2 £)- ^ e -Km)tan-M-cot^^ cot e) {m) 

The complete angular wave function now can be labeled by / and m (as a tribute to the 
spherical harmonics) and is given by 

Z?(6,ip) = ^ n =Km)-m(sinh- 1 (-cot(0)))e i ^. (61) 

It reduces to the spherical harmonics 5^ m (#, (f) for a = b = 0. In this way, the Romanovski 
polynomials shape the angular part of the wave function in the problem under consideration. 
In the following, we shall refer to Z] n (9 ) ip) as "non-spherical angular functions". 
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In Fig. 5 we display two of the lowest \ Z™ (8, <p) | functions for illustrative purposes. A more 
extended sampler can be found in Ref. 42|. A comment is in order on \Z] n (8,ip)\. In that 
regard, it is important to become aware of the fact already mentioned above that the Scarf 
II potential possesses su(l, 1) as a potential algebra, a result reported by Refs. I], [m] among 
others. There, it was pointed out that the respective Hamiltonian, H, equals H = —C — |, 
with C being the su(l, 1) Casimir operator, whose eigenva 
with j > versus + and j < in the convention of 
state solutions to Scarf II are assigned to infinite discrete unitary irreducible representations, 
{D+("O(0,p)}, of the SU(1,1) group. The SU(1, 1) labels m', and j are mapped onto ours 
via 



ues in our convention are j(j — 1) 
34j . As a consequence, the bound 



rn 



a + 



j = m — n , 



m 



j,j + l,j + 2, 



(62) 



1=1(1+1)+!, 

meaning that both j and m! are a half-integer. The representations are infinite because for 
a fixed j value, m! is bound from below to m' min = j, but it is not bound from above. 

In terms of the SU(1, 1) labels the energy rewrites as e n = — (j — |) 2 . The condition 
a > n translates now as j > \. In result, 6(0) becomes 



0(8) = ^ n=rn >-j (sinh _1 (- cot 8)) = J (1 + cot 



-2btan- 1 (-cote)Jl(' 1 ^' ~ 2b \— cot 



m'—j 



(63) 



Here we kept the parameter b general because its value does not affect the SU(1, 1) sym- 
metry. Within this context, I^W-j (sinh _1 (— cot 8) \ can be viewed as the absolute value of 
a component of a irreducible SU(1, 1) representation 40|, 41] realized in terms of the Ro- 



manovski polynomials. The \Zj n (8,tp)\ functions are then images in polar coordinate space 

components. 



ofthe| J D + (m'=i{i+i)+\) 

3=' m +2 



A. Romanovski polynomials and associated Legendre functions. 

It is quite instructive to consider the case of a vanishing 1^(9), i.e. c = 0, and compare 
eq. fl54l to eq. (j7]) for b — 0. In this case 

l = a , m 2 = (l-n) 2 , (64) 
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Fig. 5. Graphical presentation of the non-spherical-angular functions \Z\{9, tp)\ (left), and 
IZ^O, tp)\ (right) to the V^r, 9) potential in eq. (j33J). They portray in polar coordinate space in 
turn the components \D + _ 3 ( m ~z)(0,(p)\, and \D + _ a ( m ~~\o,(p)\ of the respective infinite 

discrete unitary SU(1, 1) representation. 

which allows one to relate n to I and m as m = I — n. As long as the two equations are 
equivalent, their solutions differ at most by a constant factor. This allows to establish a 
relationship between the associated Legendre functions and the Scarf II wave functions. 
Considering eqs. (124)) . and (1331) together with eqs. (1521) . one finds cot0 = — sinhz which 
produces the following intriguing relationship between the associated Legendre functions 
and the Romanovski polynomials 

P^icosd) ~ (1 + cot 2 0)-*i2, ( itr O) (- cot9 )> l-m = n = 0,l,2,...l. (65) 

In substituting the latter expression into the orthogonality integral between the associated 
Legendre functions, 

/ i?*(cos0)i^(cos0)dcos0 = 0, l^l', (66) 
Jo 

results in the following integral 

/V + cot 2 0)-^i?^'°Vcot0)4-t?' O) (-cot0)dcos0 =0, l^l'. (67) 
Jo 

Rewriting in conventional notations, the latter expression becomes 

co l -r x 

x = sinh z , I — n = I' — n = m > . (68) 
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Notion 


Symbol 


w(x) 


Interval 


Number of orth. polynomials 


Ta rrtni 


r n \X ) 


(-\ — tW1 -I- t)V 


1 — 1 11 


lnfiTiitp 

llllilll IjC 


Hermite 


H(x) 


-X 2 

e x 


[ — OO, Cxo] 


infinite 


Laguerre 


L a ' P (x) 


x P e -otx 2 


[0,oo] 


infinite 


Romanovski 




(1 + x 2 )- p e 9tan " lx 


[ — OO, CXD] 


finite 



Table 1. Characteristics of the orthogonal polynomial solutions to the generalized hypergeometric 
equation. 

This integral describes orthogonality between an infinite set of Romanovski polynomials 
with different polynomial parameters (they would define wave functions of states bound in 
different potentials) . This new orthogonality relationship does not contradict the finite or- 
thogonality in eq. (!2~lT) which is valid for states belonging to same potential (equal polynomial 
parameters). Rather, for different potentials eq. ( 1211) can be fulfilled for an infinite number 
of states. To see this let us consider for simplicity n — n' — l—m, i.e., I = I'. Given p = I + 1, 
the condition in eq. (I2T!) defines normalizability and takes the form 

2(7 -m) < 2(1 + -) - 1 = 21, (69) 

which is automatically fulfilled for any m > 0. The presence of the additional factor of 
(1 + x 2 ) -1 guarantees convergence also for m = 0. Equation (|68|) reveals that for parameters 
attached to the degree of the polynomial, an infinite number of Romanovski polynomials can 
appear orthogonal, although not precisely with respect to the weight function that defines 
their Rodrigues representation. The study presented here is similar to Ref. [43j. There, the 
exact solutions of the Schrodinger equation with the trigonometric Rosen-Morse potential 
have been expressed in terms of Romanovski polynomials (not recognized as such at that 
time) and also with parameters that depended on the degree of the polynomial. Also in this 
case, the n-dependence of the parameters, and the corresponding varying weight function 
allowed to fulfill eq. (|2~T1) for an infinitely many polynomials. 
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V. SUMMARY 



In this work we presented the classification of the orthogonal polynomial solutions to the 
generalized hypergeometric equation in the schemes of Koepf-Masjed-Jamei [221 ]. on the one 
side, and Nikiforov-Uvarov [7|], on the other. We found among them the real polynomials 
that define the solutions of the bound states within the hyperbolic Scarf potential. These 
so called Romanovski polynomials have the remarkable property that for any given set of 
parameters, only a finite number of them is orthogonal. In such a manner, the finite number 
of bound states within Scarf II were mapped onto a finite set of orthogonal polynomials of 
a new type. 

We showed that the Romanovski polynomials define also the angular part of the wave 
function of the non-central potential considered in section IV. Yet, in this case, the polynomial 
parameters resulted dependent on the polynomial degree. We identified these non-spherical 
angular solutions to the non-central potential under investigation as images in polar coordi- 
nate space of components of infinite discrete unitary £77(1, 1) representations. In the limit 
of the vanishing non-central piece of the potential, we established a non-linear relationship 
between the Romanovski polynomials and the associated Legendre functions. On the basis 
of the orthogonality integral for the latter we derived a new such integral for the former. 

The presentation contains all the details which in our understanding are indispensable 
for reproducing our results. With that we worked out two problems which could be used in 
the class on quantum mechanics and on mathematical methods in physics as well and which 
allow to practice performing with symbolic softwares. The appeal of the two examples is 
that they simultaneously relate to relevant peer research. 

The hyperbolic Scarf potential and its exact solutions are interesting mathematical enti- 
ties on their own, with several applications in physics, ranging from super-symmetric quan- 
tum mechanics over soliton physics to field theory. We expect future research to reveal more 
and interesting properties and problems related to the hyperbolic Scarf potential and its 
exact real polynomial solutions. 
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